Dynamical analysis of the buildup process near resonance 
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The time evolution of the buildup process inside a double-barrier system for off-resonance incidence 
energies is studied by considering the analytic solution of the time dependent Schrodinger equation 
with cutoff plane wave initial conditions. We show that the buildup process exhibits invariances 
under arbitrary changes on the system parameters, which can be successfully described by a simple 
and easy-to-use one-level formula. We find that the buildup of the off-resonant probability density is 
characterized by an oscillatory pattern modulated by the resonant case which governs the duration 
of the transient regime. This is evidence that off-resonant and resonant tunneling are two correlated 
processes, whose transient regime is characterized by the same transient time constant of two 
lifetimes. 

PACS: 03.65,73.40.Gk 



The buildup process of electrons inside the quantum 
well of a double barrier (DB) resonant structure has been 
one of the most important problems under investigation 
since it governs-the ultimate speed of high frequency tun- 
neling devicestlH. Although the first theoretical efforts to 
estimate the relevant time scales fcg this mechanism were 
based on stationary approachesH'B, it has been wic 
recognized that this process is of a dynamical nature^ 
hence, the solution of the time-dependent Schrodinger 
equation provides the most reliable way to tackle this 
fundamental problem. 

In this letter we provide a full quantum dynami- 
cal study of the buildup process at resonance and off- 
resonance incidence energies, based on an exact analyti- 
cal solution. We show that within the complexity of the 
process, underlying invariances can be found in the time 
evolution of the probability density. Important conclu- 
sions on the relevant time scale that characterizes both 
the resonant and off-resonant buildup are brought out 
from such invariances. 

Our analysis deals with the analytic solution of the 
time-dependent Schrodinger equation for a finite range 
potential V{x) that vanishes outside the region < x < 
L, with a cutoff plane initial condition, 
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* {x,k;t = 0) = 



— OO < X < 0, 



(1) 



0, a; > 0, 

which refers to a perfepjt reflecting shutteJi. The solution 
for the internal regioncl reads, 

*(a;,/c;t) = (l){x,k)M{0,k;t) ~ (j)*{x,k)M{0,-k]t) 



where (j}{x) stands for the stationary solution, and 0„ = 
2kun{0)un{x) / {k'^ — k'^) is given in terms of the reso- 
nant states, {w„(a;)}, and the 5— matrix poles, {fc„}, 
of the system. The index n runs over the complex 
poles fc„, distributed in the third and fourth quad- 
rants in the complex fc-plane. In the above equation 
the Moshinsky functions AI (0, q; t) are defined in terms 
of the complex error function wlz), as AI{0,q;t) = 
M{yq) — 'w{iyq)/2, where the argument is given by 
Uq = -exp{—iTr/A){Tn/2h)^/'^ [{hq/m)t^/'^], and q stands 
either for ±k or fc±„. 

The formal solution for the internal region, Eq. (||), 
allows us to study both the spatial behavior and the 
time evolution of the probability density for any inci- 
dence energy E = n^k^/2m. The main ingredients of 
Eq. (|^) are the stationary wavefunction (/){x,k), the 
resonance parameters {£'„ = e„ — iF„/2 — h^k'^/2m} 
and the corresponding resonant eigenfunctions {u„(a;)}. 
The latter can be obtained by a straightforward calcu- 
lation using the transfer matrix method adapted to the 
complex eigenvalue problenJZI. For systems with isolated 
(non-overlapping) resonances, the single-resonance term 
approximation to Eq. (^ gives an excellent description, 
except for very short times {t <^ ?i/r„) in which one 
has to consider contributions from additional resonance 
terms. In all the numerical examples presented here the 
single term approximation applies. 

As a first example, let us consider the DB struc- 
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ture (system A) with parameters: barrier heights Vi = 
V2 = 0.23 eV ^ barrier widths 61 = 62 = 5.0 nm, well 
width L0{) = 5.0 nm, and effective mass for the electron 
m = 0.067TOe. The resonance parameters for the first 
resonant state are: energy position, Si = 80.11 meV ^ 
and resonance width, Ti = 1.03 meV. In Fig. |l] (a) 
we plot l^tp, calculated from Eq. (||), as a function 
of the position x along the internal region, for specific 
times {ti\ whose increasing values are given in the figure. 
Here, the incidence energy is chosen below resonance, at 
E — 75.0 meV. Note that the off-resonant buildup oc- 
curs in such a way that is found sometimes above 
or below the asymptotic value \4>\'^. This behavior is dra- 
matically different from the monotonic growth that char- 
acterizes the special case of incidence at resonance, see 
Fig. 4 of Ref. 8. In order to show the time dependence 
of the probability density for a fixed position xq , we plot 
in part (b) |\E'(xo, fc; t)p versus t for different deviations 
from resonance AE = \E-ei\. Note that jl-p fluctuates 
around its asymptotic value |(/>p, and exhibits an oscilla- 
tory behavior not present in the resonant case, as shown 
in Fig. I (b). 

Up to here, we have illustrated the behavior of the off- 
resonance buildup only for a particular potential profile. 
It is clear that any changes in either the incident energy 
or the potential profile parameters will affect the solution 
^(x,fc;t) since the relevant input to Eq. (||), namely, 
(j){x,k), Un (x) and i?„, strongly depends on the poten- 
tial parameters. For instance, let us consider two addi- 
tional DB systems with potential profiles quite different 
from system A; the first corresponds to the symmetrical 
structure (system B) with parameters: barrier heights 
Vi ~ V2 = 0.5 eV, barrier widths 61 = 62 = 3.0 nm 
and well width loo — 10.0 nm; the second corresponds to 
an asymmetrical structure (system C), whose parameters 
are: Vi = 0.45 eV, V2 = 0.35 eV, bi = 3.0 nm, 62 = 10.0 
nm and uq = 8.0 nm. The resonance parameters for the 
first eigenstate are: Si = 37.80 meV, Fi = 0.12 meV 
(system B); and ei — 51.29 meV, Fi = 0.17 meV (sys- 
tem C). The value of the transmission peak T{ei) is unity 
for A and B; for system C, T{ei) < 1, since it is asymmet- 
ric. We choose here the incidence energies E = ei — AE 
such that the ratio 7 = T{E)/T{ei) is the same for A, 
B, and C. For example, if we choose E such that T{E) 
is 1.0 % of T{ei), from numerical inspection from a T 
versus E plot (not shown here), the incidence energies 
for A, B, and C, must be 74.97, 37.20 and 50.44 (meV), 
respectively. The results of the comparison of the time 
evolution of j^*!^ for this selection (7 — .01), are shown 
in Fig. H (a). The three curves are strongly different, 
as expected. Thus, the complete characterization of the 
buildup process for a broad range of potential geometries 
seems to be a too involved task; however, one of the pur- 
poses of this work is to show that despite this complex 
situation, the probability density has striking invariances 
under changes in the potential profiles. To illustrate the 
above let us consider a more suitable representation for 
the probability density i. e. \^ /4'\^ as a function of r, 



which is the time normalized to lifetime units. We find 
a striking result: all curves coincide exactly for the three 
different systems, see Fig. ^ (b). This result suggests the 
existence of an underlying invariance in the process. 

In order to show the existence of such invariance, we 
derive a one-level formula for the normalized probability 
density starting from the formal solution (^ . We pnpceed 
along the same lines discussed in our recent workO, but 
considering incidence energies different from resonance 
[E 7^ e„). Following such a procedure we obtain. 



I* (t) = 1 + e-^ - 2e-^/^ cos [w„t] , 



(3) 



where u;„ = (£„ — E)/Tn is a dimensionless frequency. 
The reliability of this one-level formula is shown in a plot 
of Eq. (^) included in Fig. ^ (b) , and we see an excellent 
agreement. Furthermore, note that Eq. (^ depends in 
general on the system parameters through the frequency 
LOn; however, it can be shown straightforwardly that the 
condition previously imposed on the ratio T(i?)/T(e„) 
guarantees the independence of Eq. (^ on the poten- 
tial profile. Consider the Breit-Wigner expression for the 
transmission coefficientQ, 



T{E) 



{E 
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where r° and F^ are the partial decay widths of the sys- 
tem which satisfy r„ = FO F^. This formula of T{E) 
is valid for isolated non-overlapping resonances, that is 
Tn ^ \£n — £n±i|, which is the case of a broad range 
of typical DB structures. From Eq. (P we can easily 
calculate the ratio 7 = T{E) /T{ei) and obtain an ex- 

1 /2 

pression for the frequency uji = (7"^ — l) /2. As a 
consequence, Eq. (^ is no longer dependent on the sys- 
tem parameters since it only depends on 7. Note that 
the frequency oji also measures the deviations of the in- 
cidence energy E from the resonance ei in multiple num- 
bers of Fi, i.e. AE = wiFi. In other words, different 
systems share the same curve of the probability density 
provided that deviations from resonance are the same in 
units of the corresponding resonance width Fi. Note also 
that Eq. (^) is independent of the choice ±a;„ , which im- 
plies that deviations above and below resonance give the 
same result. Since in our example 7 = 0.01, we have 
that bJi « 5.0; this can also be verified by computing the 
values of uji = |Ai?|/ri from the incidence energies used 
in Fig. |. 

Note another interesting regularity of j^* (t) /0| ; the 
damped oscillatory behavior in Eq. (|^) is modulated by 
the lower envelope 



l*(r)/0r=( 



1 



-t/tq 



(5) 



which is exactly the capacitor-like buildup law obtained 
for the special case of incidence at resonanceEl, where the 
transient time constant tq of the process is exactly two 
lifetimes, tq = 2. A plot of Eq. (||) is included in Fig. 
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^ (b). This result is relevant from a physical point of 
view, since it is a manifestation of the subtle interplay 
between the incident off-resonant carriers and the qua- 
sibound state of the system: resonant and off-resonant 
buildup, although different processes, arc not uncorrc- 
lated at all, the latter is governed by the former in the 
way exhibited in Fig. 2(b). As a consequence, the tran- 
sient regime for both situations is characterized by the 
same transient time constant tq. The above mentioned 
quantity is relevant for the design and optimization of 
resonant tunneling diodes; in this respect, a detailed dis- 
cussion caa be found in a recent work by Luryi and 
Zaslavskyll3, in which the distinction between capacitive 
and quantum contributions to an effective time constant 
is analyzed. 

In conclusion, the dynamics of the buildup mechanism 
at off-resonance incidence energies has been explored in 
typical DB resonant structures. We have shown that, 
despite the complexity that characterizes the dynamical 
process, the time evolution of the probability density ex- 
hibits invariances under arbitrary changes on the system 
parameters. From such invariances we conclude that the 
transient regime in both resonant and off-resonant pro- 
cesses is characterized by the same transient time con- 
stant of two lifetimes. Our results are valid for any DB 
system with isolated resonances and incident plane wave 
initial condition. 
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FIG. 1. (a) The birth of I'^f inside the structure as a func- 
tion of the position x, for increasing values of time: ti — 0.04 
ps, t2 = 0.4 ps ta = 0.8 ps, and = 1.2 ps (solid lines). The 
stationary solution \(fif (dashed line) is also included for com- 
parison, (b) The time evolution of I'J'f at the fixed position 
at the center of the well, for different values of AE — AEk, 
where: A^i = 0.6, A£2 = 1.1, and AE3 = 1.6 (meV). 

FIG. 2. (a) The time evolution of in the center of 
the well at off-resonance incidence energy using Eq. (2) for 
systems A, B, and C. (b) Also from Eq. (2), shows the time 
evolution of [^^(r)/^^ as a function of the time r given now in 
lifetime units; the curves of A, B and C are indistinguishable 
among them. The calculation using the one-level formula, Eq. 
(3), is also included in (b) for comparison and is indistinguish- 
able from A, B, and C. The lower envelope calculated from 
Eq. (5) is also shown. 



Also at Facultad de Ciencias, Universidad Autonoma de 
Baja California. Apartado Postal 1880, Ensenada, Baja 
California, Mexico. E-mail: villavicQfaro.ens.uabc.mx. 

1 T. C. L. G. SoUner, W. D. Goodhue, RE. Tannenwald, C. 
D Parker and D. D. Peck, Appl. Phys. Lett. 43, 588 (1983). 

^ S. Luryi, Appl. Phys. Lett. 47, 490 (1985). 

^ B. Ricco and M. Ya Azbel, Phys. Rev. B 29, 1970 (1984). 

* T. C. L. G. SoUner, E. R. Brown, W. D. Goodhue, and H. 
Q. Le, Appl. Phys. Lett. 50, 332 (1987); M. Tsuchiya, T. 
Matsusue, and H. Sakaki, Phys. Rev. Lett. 59, 2356 (1987); 
H. Yoshimura, J. N. Schulman, H. Sakaki, Phys. Rev. Lett. 
64, 2422 (1990); M. A. Talebian and W. Potz, Appl. Phys. 
Lett. 69, 1148 (1996). 

^ Note that the shutter is a device that aids to visualize the 
initial condition and hence it is not part of the system. 

^ This solution was obtained by one of the authors (J. V.) as 
an extension of the solution for the absorbing shutter in- 
troduced by G. Garcia-Calderon and A. Rubio, Phys. Rev. 
A 55, 3361 (1997). 

G. Garci'a-Calderon, R. Romo and A. Rubio, Phys. Rev. B 
50, 15142 (1994). 

* G. Garci'a-Calderon and A. Rubio, Phys. Rev. A 55, 3361 



3 



Figure 1 (a) 
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Figure 1 (b) 
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